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MICROLENSING OF A BICONICAL BROAD LINE REGION 

C. Abajas 1 , E. Mediavilla 1 , J. A. Munoz 2 , P. Gomez- Alvarez 1 and R. Gil-Merino 3 

ABSTRACT 

The influence of microlensing in the profiles of the emission lines generated in 
a biconical geometry is discussed. Microlensing amplification in this anisotropic 
model is not directly related to the bicone's intrinsic size but depends on the 
orientation of the bicone axis and on the cone aperture. The orientation of the 
projected bicone with respect to the shear of the magnification pattern can induce 
very interesting effects, like the quasi-periodic enhancements of the red/blue part 
of the emission line profile or the lack of correlation between the broad line 
region (BLR) and continuum light curves of QSOs. The emission line profiles 
of a BLR moving in a high caustic concentration exhibit sharp features that are 
well defined in wavelength. These features (spikes) correspond to the scanning 
of the kinematics of the BLR by the caustic clusters. The biconical model can 
qualitatively reproduce with a transversal (with respect to the shear) movement 
of the BLR, the recurrent blue-wing enhancement detected in the emission line 
profile of the A image of the quasar lensed system SDSS J1004+4112. The 
probability of observing this repetitive event is almost a 2% for a fraction of 
matter in stars of a 5%. This result would make plausible the detection of the 
spectral variability in SDSS J1004+4112 under the hypothesis of microlensing of 
a bicone. 
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1. INTRODUCTION 

The variability of the light-curve of a lensed object induced by the discretization of 
the mass distribution of the lens galaxy in stars or other compact objects (gravitational 
microlensing, Chang & Refsdal 1979; Young 1981; Paczynski 1986; Kayser, Refsdal & Stabell 
1986; Schneider & Weiss 1987; Irwin et al. 1989; Witt 1990; Wambsganss 2001; Gil-Merino 
& Lewis 2005) offers a promising way of studying the unresolved source of the continuum of 
QSOs. Other characteristic regions of QSOs and AGN, such as the broad line region (BLR) 
and the narrow line region (NLR), were considered too large (Rees 1984) to be significantly 
affected by the microlensing effect (Nemiroff 1988; Schneider & Wambsganss 1990). However, 
recent studies by Abajas et al. (2002) and Lewis & Ibata (2004) have shown that, according 
to the new estimates by Wandel, Peterson & Malkan (1999), the size of the BLR is small 
enough to be significantly amplified by microlensing of stellar-sized objects. Other AGN 
components such as the broad absorption line region may also be studied by microlensing 
(e.g. Belle & Lewis 2000; Chelouche 2005). 

In studies by Abajas et al. (2002) and Lewis & Ibata (2004) a number of plausible 
geometries for the BLR, taken from the complete analysis of Robinson (1995), have been 
considered. These models include, in decreasing order of symmetry, spherical, cylindrical 
and biconical geometries, which are combined with different emissivity and velocity laws. 
Additional models for the BLR could be found in Murray et al. 1995. An expected result 
confirmed by Abajas et al. (2002) and Lewis & Ibata (2004) is that the less symmetrical 
models induce the most interesting distortions in the line profiles. However, the microlensing 
of the biconical model has been studied in only two limiting cases, for inclinations 0° and 90°, 
in which the crucial effects of differential amplification between the receding and approaching 
parts of the velocity field are absent. Another important characteristic of the biconical model 
is that the length of the cones is not a direct estimate of the size of its projected surface and, 
hence, that the microlensing (de-)amplification can be strongly dependent on the orientation 
and aperture of the cones. 

In this paper we explore microlensing effects on a biconical model consistent with much 
experimental evidence (Marziani, Calvani, & Sulentic 1992; Wilson & Tsvetanov 1994) and 
theoretical work supporting the idea that the flow of emitting gas in AGN is anisotropic, 
preferentially confined to a pair of oppositely directed cones (Zheng, Binette, & Sulentic 
1990; Robinson 1995). We consider the microlensing effects at high optical depth on the 
profiles of the broad emission lines (BEL). In §3 the effects on the continuum are compared 
with those on the BLR. We then analyze the possible use of this model in explaining the 
characteristic variability detected in the BEL of the quadruple quasar lensed system SDSS 
J1004+4112. 
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2. MICROLENSING EFFECTS AT HIGH OPTICAL DEPTH ON THE 
PROFILE OF BROAD EMISSION LINES 

Throughout this paper we use the standard notation for gravitational lensing (Schneider, 
Ehlers, & Falco 1992). All microlenses have a common mass, 1M , and all distances are 
normalized to the Einstein radius associated with this microlens mass. The Einstein radius, 
i]o, projected onto the source plane is: 



where G is the gravitational constant, M is the single mass of the stars, c is the speed of 
light, and D s and D^s are the angular diameter distances to the lens, to the source, and 
from the lens to the source respectively. The Einstein radius projected onto the lens plane 
is £o = VoD d /D s . 

A standard cosmology is assumed with H = 70 kms _1 Mpc _1 , = 0.3, and il\ = 0.7. 



We have selected as an illustrative example for a general study of the microlensing of 
a biconical BLR the gravitational lens discovered by Jackson et al. (1995), B 1600+434, a 
double-imaged quasar at z s = 1.59, lensed by an edge-on spiral galaxy at z\ = 0.42. This is a 
well studied lens system that Abajas et al. (2002) found to be a good candidate for observing 
microlensing of the BLR. Given the low number of constrains for this double-imaged quasar 
we just fitted a singular isothermal sphere plus an external shear to estimate the microlensing 
parameters, the convergence or optical depth, n, and the shear due to the external mass, 7. 
For A and B images the values for («, 7) are (0.30,0.16) and (0.80,0.72) respectively. 

The normalized lens equation for a field of point masses with an external shear, 7, and 
a smooth mass distribution, k c , where the coordinate frame has been oriented such that the 
shear acts along one of the coordinate axes, is 



where x and y are the normalized image and source positions, respectively, and Xj and rrii 
are the normalized positions and masses of microlenses. The total surface mass density or 
optical depth can be written as k — + k c , where represents the contribution from the 
compact microlenses. 




(1) 



2.1. Magnification patterns 
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A ray-shooting method was used to compute the magnifications maps (Kayser, Refsdal 
& Stabell 1986; Wambsganss 1990). The map of the B1600+434B image (Fig. [TJ, with an 
area of 16?7o x I6i] at the source plane, was computed covering an area of 60£ x 60£ in the 
lens plane, and was obtained at a resolution of 1000 x 1000 pixels. 256 rays were shot per 
unlensed pixel, which means 3.6 • 10 9 rays traced between the observer and the source. The 
number of stars considered was 912, assuming that the entire optical depth is in compact 
objects, k c = 0. The theoretical mean amplification < fi t h >= |(1 — k) 2 — 7 2 | _1 is 2.068, and 
the mean amplification of the map is < /i >= 1.900. 



2.2. Biconical Broad Emission Line Region 

In this section we adopt the kinematic model for a biconical BLR used by Robinson 
(1995) to study the effects of microlensing over the line profile shapes in the context of a 
simple parameterization of some basic properties of the BLR, such as emissivity, velocity 
law, the cone aperture and inclination. 

We compute the broad emission-line profile, which depends on the wavelength A, from 
the expression 

Fx = J v e(r) 5 [x - A (l + ^) ] //(f) dV, (3) 

where \x is the amplification over the caustics pattern, 5 is the delta function, V is the total 
volume, c is the speed of light, and the emissivity and magnitude of the velocity follow radial 
dependences as: 

= (4, 

and 

v(r) = v C^~J ■ ( 5 ) 
We adopt inner (r- m ) and outer (r ex ) radii for the BLR. 

In a biconical model we need three polar coordinates (r, 9, if) measured with respect to 
the cone axis to express the projected line-of-sight velocity corresponding to an emitter, 

\ p 



v n =v Q l—) £±,|| (p>0), (6) 
with 

= sin 6* sin y> sin i + cos # cos (7a) 

£ii = sin # cos <£> sin i + cos # cos i, (7b) 
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where i is the inclination of the cones axis with respect to the line of sight, and £j_ and £y 
are the projections of the bicone when it turns around the x and y axes, respectively. Both 
projections fit when % = 0°. Figure [2] shows the perpendicular projection on the top and the 
parallel projection on the bottom. 

Following the Appendix |Al the line profile is obtained by integrating 



F x 



2ir 
I 



1+ I 

7T-0 C 



^^tt sm6d9 ) dtp, for x m |£j_,||| < \%\ < |£j 

in other cases. 



where x = c(X- X )/(v max X ), a = (e rf n c) / \X pv ) , x m = v /v v 
aperture half-angle, rj = (/3 + 3 — p)/p and 



(r in /r ex ) p , 9 C is the cone 



f=0 



f=0 



sin# cos tp, [r] ± (sin 9 sin <p cos i — cos 9 sin i)) 



[fi(r)}\ = p. {[ r }\\ °(sin^ cos^ cosi — cos 6* sin i), [r],< sin 9 sin ip) . 



(9a) 
(9b) 



2.3. Estimate of the BLR Radius 

To estimate the BLR radius for B1600+434 we apply the relationship tblr oc L 67 
from Kaspi et al. (2005) using NGC 5548 as a reference object. We have adopted the BLR 
size for high ionization lines of NGC 5548 as 2.5 light-days (He II line, Peterson & Wandel 
1999) and from NED the visible magnitude and redshift of NGC 5548, my = 13.3 and 
z = 0.017. To transform apparent magnitudes into intrinsic luminosities we have made use 
of the equation L = 47rD 1 2 um S'(l -l-^)^ -1 ^, which relates the absolute luminosity of the source, 
the luminosity distance, -Di um , the apparent flux, S, and the spectral index, a, defined by 
F v oc v~ a . We have considered a spectral index a = 0.5 (e.g., Richards et al. 2001). This is 
the procedure followed by Abajas et al. (2002) where, from an estimated intrinsic magnitude 
for B1600+434 of my = 22.69, the BLR radius for this double-imaged quasar is tblr = 4.95 
1-days for high ionization lines. 

This is a reverberation mapping size obtained, through Kaspi et al. (2005) relationship, 
from observations of many objects. In principle, we could suppose that the BLRs of these 
objects are bicones with apertures ranging from to 180 degrees (spherical case) and with 
random orientations. If the bicone axis is pointing towards the observer the reverberation 
size will depend on the bicone aperture and will approximate zero for small apertures. If 
the bicone axis were perpendicular to the observer line of sight, the reverberation sizes will 
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be luminosity-weighted radii (Koratkar & Gaskell 1991) that correspond (Robinson & Perez 
1990) in the case of a biconical geometry to 

_ J v re(r)dV _ r^Hnjl - cos O^dr _ p + 3 1 - (r in /r ex f+ 4 
rBLR I v < r ) dV J r r ™rHn{l-cos6 c )r2dr (3 + 4 1 - (r- m /r ex )P+ 3 ^ 1 j 

that is, independent of the cone aperture and, hence, coincident with the spherical case. In 
principle we should average over all the possible inclinations and cone apertures of the objects 
considered to derive the relationship of Kaspi et al. (2005) to infer the outer geometrical 
radius of the BLR (r ex ). However, having into account that the bicone orientation perpen- 
dicular to the observer is the most probable one and that we lack on information about the 
bicone aperture we finally adopt the perpendicular case (that includes the spherical case) 
with an emissivity law j3 = —1.5 and r in = 0.1 r ex . Then r ex = tblr/0.618 = 0.37i] Q , with 
l77o(lM Q ) = 21.58 1-days in this gravitational lens. This is, thus, an estimate of the intrinsic 
geometric size (probably slightly underestimated) of the BLR of one object of luminosity 
equal to that of B1600+434. Depending on the orientation and aperture of the bicone, this 
intrinsic size will result in different projected areas in the observer line of sight and, hence, 
in different microlensing amplifications. 



2.4. Microlensed Line Profiles 

Using Eq. 8 we have computed microlensed line profiles at different positions on the 
magnification pattern (Fig. [1]). Different values for the inclination, i = 0°,45°,90° and the 
cone aperture half-angle, 8 C = 10°, 30°, are considered, as well as two different exponents for 
the velocity law, p = 0.5, 2. 

Figure [3] corresponds to the i = 0° case with a cone half-aperture of 9 C = 30°. The 
exponent of the velocity law, p, is 0.5 and 2.0 for the left and right panels, respectively. In 
this limiting case the bicone is projected along its axis and would appear projected as a circle 
of radius r ex sin9 c . We obtain both strong amplification and significant deamplification of 
the line emission at different positions in the source plane. 

Figure H] corresponds to the i = 90° case with a cone half-aperture of 8 C = 30° (notice 
that now the main component of the velocity is perpendicular to the line of sight and 
the line profile has only one peak). In this other limiting case the bicone is projected 
perpendicular to its axis and would appear as a "bi-triangle" . It is evident that for this 
inclination microlensing effects on the line profile depend on the orientation of the bicone 
axis with respect to the caustics when the source moves across the magnification pattern. 
In the absence of any privileged direction in the magnification pattern this would result in a 
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random effect. However, in presence of shear, the orientation of the cone axis induces very 
interesting effects. In Fig.Hl for instance, a greater number of strong amplification events are 
found when the cone axis is parallel to the shear (lower panels) than when it is perpendicular 
(upper panels). 

In Fig. [5] the bicone has an inclination of % = 45°. For the limiting geometries considered 
before (i = 0° and % = 90°) there are no asymmetries induced by the microlensing in the line 
profile (Figs [3] and S]), assuming that obscuration and optical depth effects are negligible (see 
Murray et al. 1995). However, for any intermediate inclination microlensing induces relative 
enhancements between the blue and red parts of the profile depending on the position over 
the magnification map. It is even possible to obtain a line profile with one part of the line 
amplified and other part deamplified. Notice that microlensing effects also depend on the 
orientation of the projected cone axis with respect to the shear. 

In the previous cases we have considered a rather wide cone aperture (2 • 9 C — 60°). In 
the case of a narrower cone aperture (2 • 9 C = 20°; see Fig. [6]) the microlensing effects on the 
line profile can become extreme. 

For the cases i = 45°, 9 C = 30°, p = 0.5,2.0 and both bicone projections respect 
to the shear, we have also computed the microlensed line profiles for low ionization lines 
(corresponding to a rBLR=42.02 1-days, which means a r ex = 3.15^o)- They are shown in 
Figure [3 The effects are less important but still could be detected at some positions. 



To obtain some statistical results we have computed emission line profiles corresponding 
to the cases % = 0°,45°,90°, 9 C = 30° and both bicone projections at the position of each 
one of the pixels of the magnification map of the B1600+434B image (Fig. [I]). For each 
line profile we have estimated the magnification, the centroid displacement, an asymmetry 
parameter and the FWHM defined respectively as 



2.5. Statistics of line profiles parameters 




(11) 



(12) 



(13) 



(14) 
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where fa is the amplified line profile and gi is the non-microlensed profile multiplied by the 
mean amplification of the map < fi >. We also define the relative FWHM as (FWHM) re \ = 



The resulting accumulated probability distributions are presented in Figures [HI l9l [TU1 and 
[TTl In general, microlensing has a statistically significant impact on the line profiles generated 
by the biconical model. For instance, for the perpendicular to the shear orientation and 
intermediate inclination (45°) case, the probability of having a magnification greater than 
1 is 38%, the probability of a centroid displacement of a 0.1 is 18%, the probability of an 
asymmetry greater than 0.2 is 17%. However, changes in FWHM seem difficult to detect. 
For the case of (90°) of inclination, the effects are less significant but noticeable. As expected, 
the effects are very small in the (0°) of inclination case. 

Finally, the accumulated probabilities corresponding to the low ionization lines case 
(with % = 0°,45°,90°, 9 C = 30° and both bicone projections) indicate that the probability 
of observing flux variations in the low ionization lines is significantly smaller for the LIL 
than for the HIL (this has been observed in SDSS J1004+4112, Richards et al. 2004b, 
Gomez-Alvarez et al. 2006). However, centroid displacements and asymmetries could also 
be detected in the LIL. 



3. COMPARISON BETWEEN THE LIGHT CURVES OF THE BLR AND 

THE CONTINUUM EMISSION REGION 

We assume a disk geometry for the region emitting the continuum radiation. We suppose 
this disk to have a uniform thickness, ft < fm, and that the angle between its axis and the 
line of sight is i. 

We adopt a brightness distribution for the continuum source, 



with a power-law index p ss = 3/2 and R ss = rcontA/lO 1 /^ 65-1 ) — 1. This power- law model rep- 
resents a rough version of the exact standard model (Shalyapin et al. 2002). The continuum 
flux is computed from 



FWHM(f t )/FWHM(g t ) 




(15) 




(16) 



Recent determinations by Kochanek (2004) (see also Shalyapin et al. 2002) indicate that 
the size of the continuum source is of about a few light days. We will adopt a continuum 
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radius of 0.096^0, i-e. a radius of 6 pixels over the magnification pattern that guarantee a good 
sampling for the continuum. This value correspond to r cont = 2.07 1-days for B1600+434. 

The microlensing light-curve of the continuum (or the BLR) is found by convolving the 
brightness profile of the source with the magnification map along the path of the source over 
the pattern. In the top-left panel of Fig. [12] several paths are plotted over the magnification 
pattern of B1600+434B. For the half-aperture of the cones we adopt 6 C = 30°. In the top- 
right panel the light-curves for the BLR and continuum with % = 0° are shown. The middle 
panels correspond to i — 45° and the lower panels to i — 90°. Figure [13] is equivalent to 
Fig.H2] but with 6 C = 10°. In these figures we see that the BLR amplification may be similar 
to or even greater than that of the continuum (BLR amplifications two magnitudes greater 
than the continuum ones can be found in the perpendicular cases). For inclinations different 
from % — 0°, the degree of correlation between both the BLR and continuum light-curves 
depends on the orientation of the projected bicone axis with respect to the shear. The 
correlation is strong when the projected bicone axis is parallel to the shear direction (left 
panels in Figs. [12] [13]) but can be substantially reduced when it is perpendicular (right panels 
in Figs. [12, [13]) . To gain insight into this study (see Lewis & Ibata 2004), we have computed 
for each position (pixel) of the source in the magnification pattern pairs representing the 
BLR and continuum amplification, (/z cont , /^blr), where (from Equations 131 and [T6]) 



jF x (fi)d\ 

" fF x (<f*>)d\' (17a) 

/W fF^ t (<H>)d\ , 1 } 

where < /i > is the mean magnification of the pattern. Note that the /^BLR/Vcont ratio is 
related to changes in the equivalent width of the broad emission line. Color-scale histograms 
representing the frequency of these pairs (Lewis & Ibata 2004) for different values of incli- 
nation and cone aperture are represented in Figures HH and [15] In Fig. [14] we confirm that 
the correlation is very strong for the case i = 0°, is less marked for the parallel to the shear 
projection of the bicone (cases i = 45°, and i = 90°), and is weaker for the perpendicular 
projection (cases i = 45°, and i = 90°). The probability of finding a microlensing event in 
a strip of /zblr = /-tcont ± 0.5 decreases from 78% in the i = 0° case to 65% in the i = 90° 
perpendicular case. In the parallel case (i = 90°), the probability of finding a microlensing 
event in this strip is 71%. Similar results can be observed in Figure 10. 
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4. IS MICROLENSING AFFECTING THE BEL IN SDSS J1004+4112? 

BEL microlensing has been tentatively detected in the quadruple large-separation lensed 
quasar SDSS J1004+4112 (Inada et al. 2003) with z s = 1.732 and z x = 0.68. Blue- 
wing enhancements have been observed in the high ionization lines in component A spectra 
(Richards, Johnston & Hennawi 2004a; Gomez-Alvarez et al. 2005). A first enhancement 
was detected in two epochs in May 03 and was later undetectable in another four epochs 
from November to December 2003 (Richards et al. 2004b). A second enhancement was 
independently detected in January 04 (Gomez- Alvarez et al. 2005), March and April 2004 
(Richards, Johnston & Hennawi 2004a). In contrast, the continuum emission has been main- 
tained unaltered throughout this period. 

Microlensing of the high ionization BEL was the hypothesis originally preferred to ex- 
plain variability (Richards et al. 2004b). However this hypothesis seems difficult to reconcile 
with the recurrence of the blue-wing enhancement in the framework of the typical disk-like 
models for the BLR: instead of a new blueward enhancement, changes of the core or the red 
wing of the line were expected in so far as the caustic moves across the BLR velocity field. 
However, this behavior would not be so strange if a biconical geometry were adopted for 
the BLR. According to what we have previously seen for the case i = 45° and orientation 
of the bicone axis perpendicular to the shear, one can imagine situations in which one of 
the projected cones (blue or red) could be amplified in more than one epoch without no- 
ticeable amplification of the other. In addition, the lack of correlation between the emission 
in the line and in the continuum is also another characteristic of the biconical model. A 
singular isothermal ellipsoid (SIE) model fitted to the positions of the four images of SDSS 
J1004+4112 reveals the presence of significant shear for image A (k = 0.39, 7 = 0.64). 

The biconical model generates a two-peaked profile that can not realistically reproduce 
the emission lines in SDSS J1004+4112. However, we can consider that the bicone gives a 
contribution to the blue and red sides of a complex emission line. Another gaseous system 
will give rise to the core of the line. This two-component model has been invoked to reconcile 
the existence of an accretion disk (that also will generate two-peaked profiles) with AGN 
observations that typically show single peaked emission lines but with substructure like 
shoulders or bumps in the blue and red parts of the line profiles (see, e.g., Popovic et al. 
2004 and references therein). Thus, microlensing magnification of one of the peaks of the 
biconical model could produce an enhancement of one of the wings or an asymmetrical 
shoulder. To add more flexibility to the model we will also consider that one of the cones 
could be strongly obscured. Thus, our aim with this very simple model is mainly to test 
whether the qualitative features of the lines variability can be reproduced with microlensing 
of a bicone or not. 
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In Figs [TBI and [T71 magnification maps for SDSS J1004+4112A corresponding to two 
different fractions of matter in stars with respect to the total mass (1% and 10%) are repre- 
sented. An area of 500£ x 500£ was considered at the lens plane to compute a magnification 
pattern covering an area of 16r/ x 16r/ (1000 x 1000 pixels) at the source plane. 256 rays 
were shot per unlensed pixel, which means 2.5 • 10 11 rays traced between the observer and 
the source. The number of stars considered was 312 and 3121 for the 1% and 10% fractions 
of matter in stars, respectively. These values correspond to k c = 0.386 (1%) and n c = 0.351 
(10%). The theoretical mean amplifications, < fi t h >— 1(1 — K ) 2 — 7 2 !" 1 ; is 23.261 and 
the mean amplifications of the maps are 23.105 (1%) and 21.611 (10%), respectively. The 
most noticeable feature of both maps (Figures [16] and [TTj) is the characteristic pattern of 
alternating high and low amplification stripes induced by the strong elongation of caustics 
in the direction of the shear. 

We consider global displacements of the BLR parallel and transversal to the shear to try 
to reproduce the main features of the spectra variability in J1004+4112. The BLR radius 
for this quadruply imaged quasar is tblr = 7.57 1-days for high ionization lines (using NGC 
5548 as a reference object). In this gravitational lens lrjo(lM & ) = 16.61 1-days and, hence, 
r ex = 0.74?7o. As in the case of B1600+434 the continuum radius is 0.096?7o (6 pixels over the 
magnification pattern) that correspond to r cont = 1.59 1-days. The emissivity law is (3 — 0. 

4.1. Parallel to the Shear Paths 

In this case relative enhancements of any of the line profile peaks without noticeable 
continuum amplification are induced in a natural way by the presence of alternating stripes 
of high and low amplification. We have computed BEL profiles with r ex = 1.65%, which 
means a more conservative value for the BLR size, tblr = 17 1-days, or that the mass of 
stars in the lens are O.2M , with 1t7o(O.2M ) = 7.43 1-days. These profiles are convolved 
with the magnification pattern corresponding to a mass in stars amounting to 1% of the 
total mass (Fig. [T6"j) . following two different tracks over the pattern. As can be seen in 
Fig. [T8l the continuum remains unaltered along the entire path. However, we face here 
the same problem as Richards et al. (2004b): the enhancement event lasts several Einstein 
radii but according to the time estimates of Richards et al. (2004b) it should last about one 
tenth of an Einstein radius (for O.1M ). In principle, to reduce the time-scale of an event 
(and also the lapse between events) we could increase the number of caustics by increasing 
the fraction of mass in stars. However the caustic elongation in the shear direction is so 
strong that the time-scales of the events are not substantially reduced. The main effect of 
increasing the mass fraction in stars is to reduce the distance between consecutive stripes in 
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the magnification patterns. 

4.2. Transversal to the Shear Paths 

If the global displacement of the BLR is perpendicular to the shear, the time-scale of the 
events depends exactly on the separation between stripes. In Fig. [TH]we have represented 
the BEL profiles with r ex = li] corresponding to 20 x 20 locations regularly distributed in 
the magnification pattern corresponding to a 10% mass fraction in stars (see Fig. [TTj) . It is 
clear that a rich phenomenology of line profile variability (including recurrent or alternating 
enhancement of parts of the line profile) can be reproduced in relatively short displacements 
across the magnification pattern. Inspectioning the magnification patterns it is possible to 
find paths to roughly reproduce two consecutive events of relative enhancement of one of 
the peaks without appreciable changes in the continuum. That considered in Fig. [201 101 
instance, implies a displacement across the pattern of 0.1 Einstein radii between the third 
and the sixth line profiles from left to right; that is, ~ 2 years with v± = 700 kms^ 1 (Oguri et 
al. 2004), according to the estimates by Richards et al. (2004b), but assuming 0.5M Q stars 
in our case. Then the time between two consecutives line profiles is 0.62 years. The time 
between events of the order of 1 year can be obtained by considering that the mean mass of 
the stars is below O.5M or by supposing a factor two greater relative velocity between the 
quasar host and the lens galaxy, or by combining both possibilities. 

To study the probability of founding repetitive line variations over short time-scales 
without considerable continuum changes we have computed and statistically studied light 
curves on the magnification patterns corresponding to different fractions of matter in stars 
(2.5%, 5%, 10%, 15%, and 20%) . We consider tracks of 20 pixels corresponding to 2 years 
(with v± = 700 kms -1 ). We have extracted 1479 tracks from each magnification pattern. In 
first place, we look for continuum light curves with variability below 0.15 magnitudes. The 
probability of founding a track with this condition for the different patterns is represented 
in Figure [5TJ In second place, we identify, among these tracks, the cases in which at least 
one of the peaks follows a cycle: enhancement-fading-enhancement (if both peaks follow the 
cycle we can consider that one of them is strongly obscured) with an umbral of variability 
of at least 0.3 magnitudes (we have studied the variability of both the flux and the peak 
intensity). The probability of founding a track with both conditions, continuum variability 
below 0.15 magnitudes and repetitive line variability greater than 0.3 magnitudes, almost 
reach the 2% for the magnification pattern generated with a 5% fraction of matter in stars 
(see Fig. |2~T1) . Thus, the probability of observing this event is small but not neglectable. 

Figure [TH also illustrates other very noticeable effects that a high density of caustics 
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(originating from the high amplification and shear of J1004+4112) can induce in the line 
profiles generated by a bicone transversally oriented and transversally moving respect to the 
shear (and, hence, to the caustics). Because of the dependence of the velocity field with 
distance, each caustic will mainly affect emitters at a given velocity, enhancing a small part 
of the line profile. When more than one caustic affects the bicone, several "spikes" can be 
seen in the line profiles (see Fig. [19]). The displacement of the bicone across the caustics 
would result in the scanning of the kinematic properties of the BLR. 

Thus, the main trends of the spectral variability observed in J1004+4112 could be 
reproduced by the microlensing of a bicone moving transversally to the shear. The similarity 
of the computed phenomenon to the observed line variability can probably be improved by 
changing the proportion of matter in compact microlenses, by adjusting some of the model 
parameters (such as continuum and BLR sizes and separations) or by considering oblique 
paths. Independently of the details, the qualitative agreement between the proposed model 
and the observations in J1004+4112 supports the existence of an anisotropic gaseous system 
(susceptible of microlensing magnification) in the BLR that can give rise to substructure like 
shoulders or bumps in the line profiles. In the adopted model the inner and outer radii of the 
bicone are 9.1 and 16.6 1-days, respectively. These values are comparable to the value inferred 
from Kaspi et al. relationship for the high ionization lines in J1004+4112 (7.6 1-days). 

Statistical studies make plausible the detection of an event with features like the ones 
observed in SDSS J1004+4112 but improbable the future detection of similar events. In the 
simulations, the probability is greater for the pattern generated with a 5% fraction of matter 
in stars. 



5. CONCLUSIONS 

We have explored the effects of microlensing at high optical depth on a biconical model 
of arbitrary inclination and aperture for the BLR. We point out some results that could be 
of interest: 

1 - For symmetric geometries the size of the BLR is directly related to its projected area. 
However, a biconical BLR of any intrinsic size can have a small projected area if the cone 
aperture is small and the cone axis points toward the observer. This implies that very large 
amplifications, even greater than that of the continuum, are possible in microlensing events 
of a biconical BLR. For intermediate inclinations of the cone axis, 0° < i < 90°, microlensing 
induces relative enhancements between the red and blue parts of the line profile. 

2 - Microlensing effects (global amplification, wavelength-dependent amplification), are 
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related to the orientation of the projected bicone with respect to the shear of the mag- 
nification pattern (this could result in quasi-periodic effects when the source moves along 
the magnification pattern). In particular, the correlation between the continuum and BLR 
light-curves can be weak when the axis of the projected bicones is perpendicular to the shear. 

3 - When the bicone is oriented transversally to a caustic, the emission line profile 
presents strong enhancement of a small wavelength region (spike). The movement of the 
bicone with respect to the caustic would scan the kinematics of the BLR. When the bicone 
lies in a caustic cluster, the emission line profile can be affected by several "spikes". 

4 - We propose that the BLR can be composed by two gaseous systems. One more or 
less symmetrical that gives rise to the core of the line and other anisotropic, the bicone, that 
contributes to the blue and red sides of the line. Under this hypothesis, the main features 
related to spectral variability detected in image A of SDSS J1004+4112 (recurrent blueward 
enhancement, absence of continuum variability, and small time-scale) are qualitatively re- 
produced by microlensing of the biconical component moving perpendicular to the shear. 
The probability of an event reproducing these features is a ~ 2% in the most favorable case, 
corresponding to a fraction of mass in stars of 5%. This probability is small but makes 
plausible the detection of one event under the hypothesis of microlensing magnification of 
a bicone. However, future repetitions of the event are improbable under this hypothesis. 
Finally, it is noticeable that we can reproduce the observed variability with a bicone of inner 
and outer radii of 9.1 and 16.6 1-days, respectively. These values are comparable to the 
expected emissivity weighted radius (7.6 1-days) according to the statistical study of Kaspi 
et al. (2005). 

We thanks the anonymous referee for valuable comments and suggestions. We ac- 
knowledge the Spanish Department of Education and Science grants AYA2004-08243-C03- 
01, AYA2004-08243-C03-02 and AYA2004-08243-C03-03 for financial support. This work 
was also supported by the European Community's Sixth Framework Marie Curie Research 
Training Network Programme, Contract No. MRTN-CT-2004-505183 "ANGLES" . 



A. BLR models 



To compute the integral in Eq. [31 we define / as: 
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and then, after integrating in the r dimension and adopting x = c (A — Xo)/(v m a X Ao) and 
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Einstein Radii 

Fig. 1. — The magnification map, n(f), for the B image of B1600+434, used here to calculate 
the changes in the shape of the line profile. Each circle represents the BLR in different 
positions over the magnification map with r cx = 0.37i] . 
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Fig. 2. — Projections of the bicone, perpendicular and parallel, from left to right with a cone 
half- aperture 9 C = 30°. 
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Fig. 3. — Biconical model with i = 0°, 6 C = 30°, p = 0.5 for the left panel and p = 2 for the 
right panel. The heavy solid line is the amplified line profile and the lighter solid line is the 
line profile normalized to mean amplification in the magnification map (Fig. [I]). 
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Fig. 4. — Biconical model with i = 90°, 6 C = 30°, p = 0.5 for the left panels and p = 2 
for the right panels. The perpendicular projections of the velocity correspond to the upper 
panels and the parallel projections correspond to the lower panels. The heavy solid line is 
the amplified line profile and the lighter solid line is the line profile normalized to mean 
amplification in the magnification map (Fig. [TJ . 
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Fig. 5. — Biconical model with i = 45°, # c = 30°, p = 0.5 for the left panels and p = 2 for the 
right panels. The perpendicular projections of the velocity correspond to the upper panels, 
the parallel projections correspond to the lower panels. The heavy solid line is the amplified 
line profile and the lighter solid line is the line profile normalized to mean amplification in 
the magnification map (Fig. [1]). 
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Fig. 6. — Same as in Fig. |3J but for 9 C = 10°. 
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Fig. 7. — Same as in Fig. but for low ionization lines (r ex = 3.15?7o). 
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Fig. 8. — Accumulated probability of relative magnification (see text) in B1600+434B for 
a biconical model. Left panel: high ionization lines (r ex = 0.37r] ). Right panel: low 
ionization lines (r ex = 3.15r/o). Red, blue and green colors correspond to i — 0°, 45° and 90°, 
respectively. Solid and dotted lines correspond to the perpendicular and parallel projections 
to the shear, respectively. 
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Fig. 9. — Accumulated probability of centroid displacement (see text) in B1600+434B for 
a biconical model. Left panel: high ionization lines (r ex = 0.37r] ). Right panel: low 
ionization lines (r ex = 3.15r/o). Red, blue and green colors correspond to i — 0°, 45° and 90°, 
respectively. Solid and dotted lines correspond to the perpendicular and parallel projections 
to the shear, respectively. 
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Fig. 10. — Accumulated probability of the asymmetry parameter (see text) in B1600+434B 
for a biconical model. Left panel: high ionization lines (r ex = 0.37?7o). Right panel: low 
ionization lines (r ex = 3.15r/o). Red, blue and green colors correspond to i — 0°,45° and 90°, 
respectively. Solid and dotted lines correspond to the perpendicular and parallel projections 
to the shear, respectively. 
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Fig. 11. — Accumulated probability of relative FWHM (see text) in B1600+434B for a bicon- 
ical model. Left panel: high ionization lines (r ex = 0.37?7o). Right panel: low ionization lines 
(r ex = 3. 1577o). Red, blue and green colors correspond to i — 0°,45° and 90°, respectively. 
Solid and dotted lines correspond to the perpendicular and parallel projections to the shear, 
respectively. 
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Fig. 12. — Light-curves of the BLR and the continuum represented in solid and lighter solid 
line, respectively. Each pair of light-curves in each panel is computed across different paths 
drawn in top-left panel, correspond to the magnification map in Fig{TJ where big circles 
represent to the BLR with r ex = 0.37r] , and small circles represent to the continuum with 
^cont = 0.0967/o. The right panels correspond to the perpendicular projection of the bicone, 
and the left panels to parallel projection. Inclination increases from top to bottom, with 
0°,45° and 90°. For all cases, 6 C = 30°. 
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Fig. 14. — The correlation between the amplification of the BLR and the continuum is 
represented, when a biconical model is considered for the BLR and a Keplerian disk for the 
continuum. The right panels correspond to the perpendicular projection of the bicone and 
the left panels to parallel projection. Inclination increases from top to bottom, with 0°,45° 
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Fig. 16. — The magnification map for the A image of SDSS J1004+4112 for 1% of total 
mass in stars. Each point represents to BLR in different positions over the magnification 
map, and the lines are the paths over the map where the light-curves are computed. The 
biggest circle represents to the BLR with r ex = 1.65??o, and the smallest circle represents to 
the continuum with r cont = 0. 096^0- 
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Fig. 17.— The magnification map for the A image of SDSS J1004+4112 for 10% of total 
mass in stars. The biggest circle represents to the BLR with r cx = I770, and the smallest 
circle represents to the continuum with r cont = 0.096i] . The small line starting from the 
circles centre is the path over the magnification pattern where the light-curves are computed. 
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Fig. 18. — The line profiles for the BLR for each point and the light-curves for the BLR 
and the continuum for each path marked in Fig. [IS] are shown. The parameter are i = 45°, 
9 C = 5°, p = 2, r cont = 0.096?7o, r ex = l-65?7o and r in = 0.lr] . In the line profile panels the 
heavy solid lines represent the amplified line profile and the lighter solid lines correspond 
to the line profile normalized to mean amplification in the magnification map. In the light- 
curve panels the solid line represents the light-curve from the BLR and the lighter solid line 
corresponds to the light-curve from the continuum. 
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Fig. 19. — The BEL profiles for 20 x 20 locations regularly distributed in the magnification 
pattern (Fig. [TTj) are shown. The parameter are % = 45°, 6 C = 5°, p = 2, r cx = I770 and 
r in = 0.1?7o. The heavy solid lines represent the amplified line profile and the lighter solid 
lines correspond to the line profile normalized to mean amplification in the magnification 
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Fig. 20. — The line profiles for the BLR and the light-curves for the BLR and the continuum 
for a track superimposed on the magnification pattern (Fig. [T71) are shown. The parameter 
are i = 45°, 6 C = 5°, p = 2 , r cont = 0.096r] , r ex = lr} and r in = 0.55?7o. In the line 
profile panel the heavy solid lines represent the amplified line profile and the lighter solid 
lines correspond to the line profile normalized to mean amplification in the magnification 
map. In the light-curve panel the solid line represents the light curve from the BLR and the 
lighter solid line corresponds to the light-curve from the continuum. 
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Fig. 21. — The probability of founding repetitive line variations over short time-scales with- 
out considerable continuum's changes in our simulated magnification maps is represented 
versus the fraction of matter in stars. The dot line represents the probability of founding a 
track where the variability of the continuum light curves are below 0.15 magnitudes. The 
dash and solid lines represent the probability of founding a track where the variability of 
the continuum light curves are below 0.15 magnitudes, and repetitive line variability greater 
than 0.3 magnitudes, taking into account the peak and the flux intensity respectively (See 
text). 



